function d = hybridNeuron_E(t, y, params)
    % hybridNeuronE: 定义含电场 E 影响的神经元模型的微分方程
    % dz/dτ 项中增加了电场影响 E，以研究外部电场对神经元动力学的作用。
    % 状态向量 y = [x; y; z]
    % params 为参数结构体，包含字段：
    %        a, b, c, alpha, beta, delta, A, f, E
    %
    % 对应论文公式(6):

    x     = y(1);
    y_var = y(2);
    z     = y(3);
    
    % 外部激励信号 u_s = A*cos(f*t)
    u_s = params.A * cos(params.f * t);
    
    % 计算各个微分方程分量
    dxdt = (1 - params.a)*x - (1/3)*x^3 - y_var + u_s;
    dydt = params.c * (x - params.b * sin(params.beta * z) * y_var);
    % 在原始模型基础上，dz/dτ 增加电场 E 的影响
    dzdt = params.delta * y_var - params.alpha * z + params.E;
    
    % 返回状态导数
    d = [dxdt; dydt; dzdt];
end
